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Despite several eflForts the nonlinear Meissner effect in d-wave superconductors, as has been discussed 
by Yip and Sauls in 1992, has not been verified experimentally in high-Tc superconductors at present. 
Here, we reinvestigate the nonlinear response expected in a d-wave superconductor. While the linear 
\H\ field dependence of the penetration depth, predicted by Yip and Sauls, is restricted by the lower 
critical field and can be masked by nonlocal effects, we argue that the upturn of the nonlinear 
coefficient of the quadratic field dependence is more stable and remains observable over a broader 
range of parameters. We investigate this by studying the infiuence of nonmagnetic impurities on the 
nonlinear response. We discuss the difficulties of observing this intrinsic d-wave signature in present 
day high-Tc films and single crystals. 



PACS: 74.25.Nf, 74.20. Fg, 74.72.-h 
I. INTRODUCTION 

In a dj.2_y2-wave superconductor, quasi-particles near 
the nodes of the gap give rise to an intrinsic nonlin- 
ear electrodynamic response. Yip and Saulscl discussed 
this and suggested that at sufficiently low temperatures 
this nonlincarity would lead to an increase in the pen- 
etration depth A which would vary as the magnitude 
of the magnetic field \H\ with a coefficient that de- 
pended upon the orientation of the field relative to the 
nodes. rnAlthough there have been various experimental 
studiesaB of the magnetic field dependence of the pene- 
tration depth, there are at present no observations which 
are in agreement with the size, field, and temperature d&i 
pendence predicted in Ref. ^. It is known that impuritiestl 
can wash out the linear H dependence, leading to an 
dependence for A and recently Li et aLLl have argued 
that nonlocal effects, except for special orientations of 
the field, will lead to a quadratic dependence below a 
cross-over field of order Hd. Thus, the magnitude \H\ 
signature of a dx2_y2-wave state in the penetration depth 
may be difficult at best to observe. ^ 

However, as discussed by Xu, Yip and SaulsH, at higher 
temperatures the temperature dependence of the coeffi- 
cient of the quadratic field term in the penetration depth 
shows a clear deviation from the exponentially decaying 
temperature behavior which would be observed in a fully 
gapped s-wave superconductor. In particular for a dx2^y2 
gap, this coefficient exhibits a T^^ temperature depen- 
dence so that contrary to the s-wave case, the strength 
of the contribution increases as T is decreased until 
it saturates due to impurities, nonlocal eflFects, or possi- 
bly, for the magnetic field oriented along the nodes, the 
crossover to the Yip-Sauls \H\ dependence. It is this as- 
pect of the nonlinear response that-. we wish to explore. 

As we have previously discussedlj, the quadratic term 



in the penetration depth leads to a nonlinear inductance 
which gives rise to third order intermodulation effects. 
This nonlinearity in fact represents a problem for su- 
perconducting communication filters and care must be 
taken to reduce its presence. However, from the point 
of view of studying the nonlinear superconducting re- 
sponse, micaiwave intermodulation effects provide a sen- 
sitive probdj^EI. Here we present results which illustrate 
the type of nonlinear dependence that can be expected. 

In Section || we study the field and temperature de- 
pendence of the nonlinear superfluid density for a clean 
d-wave superconductor. The crossover between the \H\ 
linear Yip-Sauls regime and the quadratic regime with its 
divergent coefficient are discussed. Section III illustrates 
the corresponding behavior expected to be seen in har- 
monic generation and intermodulation. In Section |^ we 
study the influence of nonmagnetic impurities in order to 
see how stable these signatures of the d-wave state are. 
We will see that the increase of the nonlinear coefficient 
at low temperatures appears to be more stable than the 
\H\ linear behavior due to its restriction by the lower 
critical field Hd . We discuss the difficulties in observing 
these signatures in present day high-Tc systems. Section 
^ contains our conclusions. 



II. THE NONLINEAR SUPERFLUID DENSITY 

The current density in a superconductor is given by 
the sum of the superfluid flow js and the quasi-particle 
backflow contribution 
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Here, /(e) — 1/[1 + exp(e/T)] is the Fermi function and 
we have taken a simple circular Fermi surface with a 
dx-^_y2 gap. As has been discussed by Li et alE nonlocal 
effects will become important for current flow along the 
antinodal direction for a magnetic field oriented perpen- 
dicular to the Cu02 planes. However, nonlocal effects are 
negligible for current flow along a nodal direction. Here 
and in the following we will therefore study this second 
case using the d-wavc angular dependence 



A (9) = A (T) sin 26. 
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Here, for convenience, we have chosen our coordinates 
such that = corresponds to the nodal direction. We 
have checked that qualitatively very similar results are 
found for current flow along the antinodal direction, if 
nonlocal effects are neglected. The superfluid density 
is defined by 
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so that 
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with jc = neAo/mvF the pair breaking current density. 
In the limit js ^ 0, we have 
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we have for the nonlinear contribution to the superfluid 
density 
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At small current flow this is related to the nonlinear 
change in penetration depth AA via 



AA ^ 1 Sus 



(9) 



It is this nonlinear part of the superfluid density that we 
will study. 



There are some simple limiting cases. As discussed by 
Yip and Sauls, when js/jc T'/Aq, then one can expand 
Eq. (g) in powers of js/jc so that 
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with 
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At low temperatures this expression yields 
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the 1/T divergence mentioned above. Alternatively, in 
the low temperature limit where 



'- = 1 - 2 In 2 — 

n \Ao 



and js/jc ^ T/ Aq we have the Yip-Sauls resultE. 
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For the calculations presented in the following we choose 
Aq /Tc = 3, a value thM fits low temperature penetration 
depth data on YBCO0. 
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FIG. 1. The nonlinear part of the superfluid density Sus/n 
as a function of superfluid current density js/jc for temper- 
ature T/Tc = 0.005 (solid line). The dashed line shows the 
parabohc low-js expansion Eq. ( |lo| ) and the dotted line the 
Yip-Sauls result Eq. (0). 
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FIG. 2. The nonlinear part of the superfluid density Sus/n 
as a function of superfluid current density js/jc for different 
reduced temperatures T/Tc. 



For T/Tc = 0.005 the two limits in Eqs. (|T^) and 
( |r^ ) are shown in Fig. Q along with the numerical resuh 
obtained by numericaUy integrating Eq. (^. Here we 
see a crossover from a quadratic dependence on js/jc for 
js/jc < T/Tc to a linear dependence when js/jc becomes 
larger than T/Tc- Fig. || shows a sequence of curves for 
Sus/n versus js/jc for different reduced temperatures. At 
larger values of js/jc, 6ns/n has approximately the same 
slope, but the curves are shifted down by an amount pro- 
portiopal to T/Tc. Af has been pointed out by Yip and 
SaulsEl and Li et al.El this hnear js/jc dependence can 
at best only be observable at low temperatures due to 
the fact that a type II superconductor will enter the vor- 
tex state, if the current density level js/jc reaches the 
lower critical field Hd/Hc — 0.01 for high- Tc supercon- 
ductors. At small values js/jc < T/Tc, ^Us/n enters the 
parabolic regime and the curvature (3 (T) given in Eg. 
(O) diverges like 1 /T upon lowering the temperatureSfl. 
At the same time the convergence radius of the Taylor 
expansion decreases like T/Tc, showing the crossover to 
the nonanalytic Yip-Sauls result at T = 0. 



III. HARMONIC GENERATION AND 
INTERMODULATION 

The nonlinear response to an applied microwave field 
such that the microwave frequency is small compared to 
the quasi-particle relaxation time can be determined from 
Eqs. d) and (|). If 
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Now, only odd frequency terms arise since n<,(js) = 
ns{~js)- For example, the third harmonic jasinSwt has 
an amplitude 
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In the 'high temperature limit', T/Tc ^ js/jc we find 
using Eq. (p^) 



J3 = J/3(r)f 
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and in the 'low temperature limit', T/Tc ^ js/jc using 
Eq. (|lj) we have 
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FIG. 3. The amplitude ja/jc of the third harmonic gen- 
erated by the nonlinear superfluid density as a function of 
the amplitude of the fundamental current density jso/jc for 
different reduced temperatures T/T^. The two limiting cases 
Eqs. (|l|) (for T/T^ = 0.001) and ^ are shown as the solid 
and the dotted line, respectively. 



Fig. ^ shows a double logarithmic plot of js/jc ver- 
sus jso/jc for three different values of T/Tc and the 
crossover from Eq. (p3) (solid line) to Eq. ( p^ ) (dot- 
ted line) is clearly seenH! Note, that due to the prefactor 
/3(T) in Eq. ( p^ ) the nonlinear response js/jc in the low 
jso/jc regime increases, when the temperature is low- 
ered. Fig. ^ shows the temperature dependence of js/jc 
for jso/jc = 0.01 along with the high temperature limit 
from Eq. (^8|). The third harmonic amplitude js follows 
the 1/T divergence until T/Tc falls below jso/jc = 0.01. 
Below that point js saturates due to the crossover to the 
Yip-Sauls limit Eq. (19). This low temperature peak 
in the third harmonic amplitude is a consequence of the 
nodes of the d-wave state and does not exist for an s-wave 
superconductorO. 
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FIG. 4. Temperature dependence of js/jc for jso/ jc ~ 0.01 
(solid line). The dashed line shows the high temperature 
limit Eq. Below T/T^ = jso/jc = 0.01, the third har- 

monic amplitude js/jc saturates due to the crossover to the 
Yip- Sauls regime Eq. (p^. 



One can also explore two-tone intermodulation which 
is important in communication filter applications. Here 



js (t) = jsi sincjii + js2 sin 0^2* 
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with LUi and uj2 close in frequency so that the intermod- 
ulation frequency 2iOi — 102 lies within the pass band of 
the filter. In this case we find for the amplitude j2uji-rcn 
of the intermodulation response at high temperaturesEiJ 
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while at low temperatures, when js2 ^ jsi 
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and when jsi < js2 
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Thus, at fixed jsi and js2 the temperature dependence of 
J2wi-cj2 '^ill have the same qualitative behavior as jslT) 
in Fig. |. 

This analysis shows that measurements of third har- 
monics and intermodulation allow to directly access 
the temperature dependence of the nonlinear coefficient 
/?(T), utilizing Eqs. (|l8|) or (pT|). In principle, this allows 
to probe experimentally whether this low temperature 
peak in P exists or not. 



IV. INFLUENCE OF IMPURITIES 

Now we wish to consider the infiuence of nonmagnetic 
impurities on the nonlinear response discussed in the pre- 
vious sections in order to see how stable the d-wave sig- 
natures in the nonlinear response are. 



In the presence of impurities the total current density 
is given bycl 
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where N±{Q, uj) is the density of states for the comoving 
and countermoving quasiparticles, respectively: 
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Here, uj^w) is the renormalized frequency aud .has to be 
determined by the sclfconsistent equationst3ll3 
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Here, c is the cotangent of the scattering phase shift and 
r the scattering rate. Using Eqs. ( p^ - (^ along with 
Eqs. (^ and (0) we can extract i5ns(j7)/n in the presence 
of nonmagnetic impurities. 
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FIG. 5. Sua/n as a function oi js/jc for different reduced 
temperatures T/Tc in the presence of nonmagnetic impuri- 
ties in the unitarity limit c = 0. The scattering rate is 
F = O.OOSTc 



In Fig. ^we show Sns{js)/n in the unitarity limit c = 
with r = O.OOSTc, a typical value in a range that has 
been used to fit low tempecature penetration depthllJ and 
thermal conductivity datallJ. Due to the finite impurity 
scattering Sng/n now stays quadratic at low js/jc down 
to r = 0. Now, the linear js/jc regime is only entered 
at higher values of js/jc ^ T/Tc- The restriction due 
to the lower critical field js/jc < Hci/Hc — 0.01 can 
make the linear js/jc regime unobservable in this case, 
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as can nonlocai effects for current flow along the anti- 
nodal directionEI. 
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FIG. 6. Temperature dependence of tfie nonlinear coeffi- 
cient I3{T /Tc) for different scattering rates V in tlie unitarity 
limit. 
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FIG. 7. The density of states N{u) in the presence of im- 
purities in the unitarity limit for F = O.OOGTc. Close to the 
Fermi level oj/Aq < 0.05 impurity states are generated. How- 
ever, over a broad region N{u)) still varies linearly with lj as 
in the clean d-wave state. 



In order to find the nonlinear coefficient /3{T/Tc) it is 
advantageous to perform the calculation on the imagi- 
nary frequency axis. Then we find 

.r;^ A^(e)Re|^^i^^^|. (28) 

Here, (!;„ are the renormalized Matsubara frequencies 
which have to be determined selfconsistently by the imag- 
inary axis counterparts of Eqs. (|2^) and (p7|): 

<D„ = u)n + r o^°_2 (29) 
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with LUn — {2n+ 1)ttT being the unrenormalized Matsub- 
ara frequencies. 

Fig. ^ shows (]{T/Tc) for different scattering rates. 
Now, the l/T divergence is cut off at low temperatures 
by the impurity scattering. However, a peak, signify- 
ing the underlying d-wave nature of the superconducting 
gap, still remains unless T/Tc becomes of the order of 
0.1. Qualitatively similar results are found for nonuni- 
tary scattering c ^ 0. The influence of impurities be- 
comes even less pronounced, though. 



As an illustration we show the density of states N{lu) 
in the presence of impurities in Fig. ^. Here one 
can see the impurity states generated at low energies 
uj /Aq < O.OSHEJ. However, at higher energies there is 
still a broad region where N{uj) varies linearly with lu. 
At not too low temperatures the nonlinear response still 
picks up this linear variation resulting in an upturn of 
the nonlinear coefficient (3{T /Tc) upon lowering the tem- 
perature. 

This signature of the d-wave state should remain ob- 
servable over a much broader range of parameters than 
the |7J|-linear regime discussed by Yip and Sauls. While 
the |i?|-linear regime is limited by the lower critical field 
T/Tc < Hci/Hc ~ 0.01, this peak in l3{T/Tc) wiU re- 
main until T/Tc becomes of the order of 0.1, making it 
a better candidate for a search of d-wave behavior in the 
nonlinear response. From the study in Ref. |^ we expect 
that this upturn of f3 should also remain visible, if one 
includes nonlocal effects for the case of current flow along 
the antinodal direction. 

Some remarks concerning the difficulties of observing 
this behavior should be made at this point. In order to 
extract the coefficient (3 from measurements of the H- 
field dependence of the penetration depth A one would 
need a very high resolution, as can be estimated from 
Eqs. (||) and (p^): even if js/jc — 0.01 and for /? we 
take a typical value of 4, we find AA/A ~ 2 x 10^"* which 
challenges existing techniquesu. A more direct way to 
measure f3 would be harmonic generation or intermodu- 
lation utilizing Eqs. ( |l8| ) and (|2l| ) as has been discussed 
in the previous section. A study of the temperature de- 
pendence of intermodulation in high- Tc (TBCCO) films 
has been done in Ref. |l^. In that study no increase in 
intermodulation was found down to the lowest tempera- 
tures of 25 K, which might not be small enough, however. 
A similar result was found in measurements of third har- 
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monic generation in YBCO film3l3. These studies also 
showed that the absolute magnitude of the nonlinear re- 
sponse is higher than expected from the intrinsic d-wave 
response discussed here. This and the study in Ref. ^ 
indicates that other sources of nonlinear behavior, as for 
example weakJink grain boundaries which act as Joseph- 
son JunctionsEa, might dominate the nonlinear response 
in present day high- Tc films. Thus, the difficulties of 
observing d-wave behavior in the nonlinear response are 
mainly related to the existence of extrinsic effects even 
in the best presently available systems. Nevertheless, we 
suggest that a temperature dependent measurement of 
harmonic generation or intermodulation in the highest 
quality single crystals available today provides the best 
hope of observing d-wave behavior in the nonlinear re- 
sponse. 



V. CONCLUSIONS 

We studied the temperature and field dependence of 
the nonlinear electrodynamic response in a d-wave su- 
perconductor. The signatures of the d-wave state are the 
|/?|-linear regime, as discussed by Yip and Sauls, and 
an upturn of the nonlinear coefficient /? in the quadratic 
regime at low temperatures. This coefficient can be di- 
rectly measured by harmonic generation and intermod- 
ulation. While the |_ff|-linear regime is limited by the 
lower critical field and can be masked by impurity scat- 
tering and nonlocal effects, the upturn of the coefficient 
f3 appears to be more stable and should remain observ- 
able over a broader range of parameters. We showed this 
explicitly by studying the influence of nonmagnetic im- 
purity scattering. It is possible, however, that in present 
day high- Tc films and single crystals extrinsic effects still 
dominate the nonlinear response, masking this intrinsic 
signature of the d-wave symmetry. 
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